Prediction of the collapse point of overloaded materials by monitoring energy 

emissions 
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A bundle of many fibers with stochastically distributed breaking thresholds is considered as a 
model of composite materials. The fibers are assumed to share the load equally, and to obey 
Hookean elasticity up to the breaking point. The bundle is slightly overloaded, which leads to 
complete failure. We study the properties of emission bursts in which an amount of energy E is 
released. The analysis shows that the size of the energy bursts has a minimum when the system is 
half-way from the collapse point. 

PACS numbers: 02.50.-r 



I. INTRODUCTION 

During the failure process in composite materials un- 
der external load, bursts (avalanches) of different mag- 
nitudes are produced, where a burst consists of simulta- 
neous rupture of several elements. At each failure, the 
sudden internal stress redistribution in the material is 
accompanied by a rapid release of mechanical energy. 
Therefore, with each burst there will be a corresponding 
energy emission burst. A useful experimental technique 
to monitor [0] such energy bursts is to measure the acous- 
tic emissions, the elastically radiated waves produced in 
the bursts P,[3|. 

Fiber bundle models, with statistically distributed 
thresholds for the breakdown of the individual fibers, are 
interesting models of failure processes in materials. They 
are characterised by clear-cut rules for how stress caused 
by a failed element is redistributed on undamaged fibers. 
These models have been much studied since they can be 
analysed to an extent that is not possible for more com- 
plex materials (For reviews, see The statistical 
distribution of the magnitude of avalanches in fiber bun- 
dles is well studied |9l4ll| , and the failure dynamics under 
constant load has been formulated through recursion re- 
lations which in turn explore the phase transitions and 
associated critical behavior in these models fl2| . 

In this article we show that the catastrophic collapse 
point of an overloaded bundle can be predicted by mon- 
itoring the energy emission rate. 

We consider a bundle consisting of a large number N 
of fibers, clamped at both ends. We study equal-load- 
sharing models, in which the load previously carried by a 
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[0] The Scottish Forth Road Bridge is supported by two main cables, 
each with 11618 fibers (wires), some of which have failed. Acous- 
tic monitoring was installed in 2006 to detect further snapping 
of the wires. (See Forth Road Bridge in Wikipedia.) 



failed fiber is shared equally by all the remaining undam- 
aged fibers [l3l - [l6j . The fibers obey Hooke's law, such 
that the energy stored at elongation x equals x 2 /2, where 
we for simplicity have set the elasticity constant equal to 
unity. Each fiber i is associated with a breakdown thresh- 
old Xi for its elongation. When its length exceeds Xj, the 
fiber breaks immediately, and does not contribute to the 
bundle strength thereafter. The individual thresholds Xi 
are assumed to be independent random variables with 
the same cumulative distribution function P(x) and a 
corresponding density function p(x). If an external load 
F is applied to a fiber bundle, the resul ting breakdown 
events can be seen as a sequential process [10( . In the first 
step all fibers that cannot withstand the applied stress 
(F/N) break. Then the stress is redistributed on the 
surviving fibers, which compels further fibers to fail, etc. 
This iterative process continues until all fibers fail, or an 
equilibrium situation with a nonzero bundle strength is 
reached. Since the number of fibers are finite, the total 
number of steps, f/, in this sequential process is finite. 
At the stress (or elongation) a per surviving fiber the 
total force on the bundle is a times the number of intact 
fibers. The expected, or average, force at this stage is 
therefore 



F(a) = N a [l-P(cr)]. 



(1) 



The maximum F c of F(o~) corresponds to the value o~ c 
for which dF/da vanishes. Thus 



1 - P(a c ) - <r c p{a c ) = 0, 
where the critical stress o~ c is defined as 
a c = F c /N. 



(2) 



(3) 



When the applied load is more than F c (or a > o~ c ), we 
call that the bundle is overloaded. 

We can study the stepwise failure process in the fiber 
bundle when a fixed external load F is applied. Then 
the initial external stress is a = F/N . Let N t be the 
number of undamaged fibers at step no. t, with No = N. 
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We want to determine how N t decreases until the failure 
process stops. With N t intact fibers, an expected number 



[N P{Na/N t )} 



(4) 



of fibers will have thresholds that cannot withstand the 
load, and consequently these fibers break at once. Here 
[X] denotes the largest integer not exceeding X. The 
number of intact fibers in the next step is therefore 



N t+ i=N-[N P(Na/N t )}, 
or, for all practical purposes, 

n t+ i = 1 - P(a/n t ). 



(5) 



(0) 



Here n t denotes the fraction of undamaged fibers at step 



t, 



n t = N t /N 



(7) 



In each burst a certain amount of elastic energy is re- 
leased [l*7 . and we consider now the energy emission 
process. 



II. ENERGY RELEASE MINIMUM 

For a burst in which the number of intact fibers is 
reduced from N t _i to Nt, all fibers with thresholds x 
between the values F/N t -i = a/n t -i and F/N t = cr/n t 
break. Since there are Np(x) dx fibers with thresholds 
in (x,x + dx), the energy emitted in this burst is given 
by 



Ei 



a/n t 



\x 2 Np{x) dx. 



(8) 



We consider external stresses that are slightly above 
the critical value, 



(9) 



where e is small and positive. Large stresses are less in- 
teresting, since the system will then break down quickly 
(See the Appendix). Simulations for a uniform thresh- 
old distribution, p(x) — 1 for < x < 1, show that 
the energy emission has a minimum at some value i_E(e), 
and that for varying loads the minima all occur at value 
close to | when plotted as function of the scaled variable 
t/tf (Fig. 1). Here tf is the number of iterations corre- 
sponding to complete bundle failure. We will 'show this 
analytically, and in addition demonstrate that the result 
is not limited to the uniform distribution. Fig. [1] shows 
that for large N and near the minimum the energy emis- 
sion E t appears to be an almost continuous function of t. 
The minimum is therefore located where the derivative 
of E t with respect to t vanishes. From © we obtain 

dE t jdt = |iV [-a 2 n^ 4 p(a/n t ) n t + cr 2 ™^ p(a/n t -i) h t - 

(10) 
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Figure 1: Energy emission Et vs. step t (upper plot) and vs. 
the rescaled step variable t/tf (lower plot) for the uniform 
threshold distribution for a bundle of N = 10 7 fibers. Differ- 
ent symbols are used for different excess stress levels a — a c : 
0.001 (circles), 0.003 (triangles), 0.007 (squares). 



where the dot denotes differentiation with respect to 
time. 

To eliminate quantities at time step t — 1 we use the 
connection which gives 



n t = ht-i p(a/n t -i) an t \. 



Insertion into (fTU|) yields 



dE 
~dt 



N 



a 3 p(a/n t ) 



H-i 



We also need nt-i- Eq. (|6|) gives 

a/nt-i =P- 1 (l~n t ), 



(11) 



(12) 



(13) 



where P 1 denotes the inverse function to P. In conclu- 
sion, 



dE 
~dt 



. N 
nt — 
2cr 



a 3 p{a/n t ) 



+ (p- 1 (l-n t )) 



(14) 
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Since fit is always negative, the minimum of energy emis- 
sion occurs when the relative number of intact fibers sat- 
isfies 

P- 1 (l-n t ) = ^a3p(a/n t )nt\ (15) 

i.e. 

1 - m = P (vW^MK -4 ) • (16) 
Now we turn to specific cases. 

III. EXPLICIT RESULTS 

We start with the simplest case, the uniform distribu- 
tion, on which the simulations in Fig. 1 were based. 

A. Uniform distribution 

For the uniform threshold distribution, P(x) = x for 
< x < 1, the condition (p~6|) takes the form 

1 - nt = \/<r 3 n~ t \ (17) 

or 

(1-ntK = ff3/2 - ( 18 ) 

For this case the maximum of the force, Eq. (jTJ) , is F = 
N/4, corresponding to the critical stress <j c = 4. For a 
small excess stress e = cr — a c = a — i, Eq. ([T5]l takes the 
form 

(l-n t )n^(i + £ ) 3/2 ^i(l + 6 e + 0( e 2 )). (19) 

To lowest order n t = ^, and one shows easily that 

n t = \ + 3e (20) 

satisfies (IT91 to first order in e. 

To find ig, the value of t corresponding to (|2"0"j) we use 
the previously derived solution of the iteration (O 

n t = \-Jl taa.(At - B), (21) 

where A = tan" 1 (2 v /e) ~ 2^/1 and B = tan" 1 (l/2 v /e). 

We see that (f2TJ)) requires tan(j4fg — B) = — 3y/e, i.e. 
tE—B/A = —3y/e/A. To lowest order we therefore obtain 

t E =B/A-\. (22) 

Moreover we see from (|2"Tj) that at for 

t = t f = 2B/A (23) 
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Figure 2: The energy emission E t vs. the scaled step vari- 
able t/tf for a bundle of N = 10 7 fibers obeying the Weibull 
threshold distribution Q26p . Different symbols are used for 
different excess stress levels a— a c : 0.001 (circles), 0.003 (tri- 
angles), 0.007 (squares). 

we have 

"*/ = 5 - tan ( B ) = 0. ( 24 ) 

signifying complete collapse of the fiber bundle. By com- 
parison between ([22]) and ([23]) we obtain 

t E = \t f -\. (25) 

Thus the minimum energy emission occurs almost 
halfway to complete bundle collapse, in agreement with 
the simulations shown in Fig. [1] 



B. Weibull distribution 

To illustrate the generality of the connection between 
the energy emission minimum and the bundle collapse, 
we now turn to a completely different threshold distribu- 
tion, viz. a Weibull distribution of index 5, 

P{x) = l-e- x \ (26) 

The critical stress for this distribution is a c = (5e) -1 / 5 = 
0.5933994. Simulations reveal a qualitatively similar 
behavior as for the uniform distribution (Fig. [5] ) . 

The condition (|16[) for minimal energy emission takes 
in this case the form 

Inrn = 5 5 / 2 cr 35 / 2 V 20 exp(~5(T 5 /2n t ). (27) 

It is straightforward to verify that for a — a c — (5e) -1 / 5 
(|27p has the solution n t — e -1 / 5 . We are interested in a 
slightly overloaded bundle, 

<7 = <7 c + e= (5e)- 1/5 + e, (28) 
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with e small and positive, and we seek the corresponding 
value of n t . Putting (|2"5)l and n t — e" 1 / 5 + 8 into the 
condition f|27f>. and expanding in e and <5, we obtain to 
first order in the small quantities 



jV5 S = 3(5 e )i/i 



Thus 5=| 5 1 / 5 e. Hence 



nt 



e-V5 + 2 5 i/5 £ 



(29) 



(30) 



is the relative number of undamaged fibers when the en- 
ergy emission is minimal. 

To find t E , the value of t corresponding to (f30|) . we 
take advantage of the ground work already done in [l8[ 
for the Weibull distribution. Eq. (29) in [18( shows that 
for small e the iteration is of the form 



= e - 1/5 - 5 1/5 v^7c* taxi(tVeC - c), 



(31) 



with C = l^e) 1 / 5 and the constant c = tan 1 [(1 — 
e" 1 / 5 ) 5 _1 / 5 \JC/e\ ensures that the initial condition 
no = 1 is satisfied. 

Comparison between (|3Tfl) and (|3"Tj) gives 

tan(t£;\/ea-c) = -^VeC. (32) 

To dominating order, then, 

t E VdJ-c= -tan-^fVd?) =s (33) 

For small e the constant c is very close to 7r/2, so that in 
good approximation we have 



t E 



2VeC* 



(34) 



The collapse time tt was also evaluated in [18[ to be 



Consequently we have 



tE 



2 l f 



(36) 



To excellent approximation the minimum of energy 
emission occur halfway to the final collapse also for the 
Weibull distribution. There is every reason to believe 
that this feature is general. Note that since t E and tf 
are large numbers, the constants | and | in (|2"5)l and 
(|5rJ|) are of no significance. 



IV. CONCLUDING REMARKS 

In summary, we have considered energy emission burst 
from slightly overloaded fiber bundles. During the degra- 
dation process there is a stage tE at which the energy 
emission is minimal, and we have demonstrated that the 
total bundle collapse occurs near 2t E - The demonstra- 
tion has been performed merely for two very different 
distributions of the fiber breaking thresholds, but the re- 
sult is doubtlessly universal. Thus the minimal energy 
emission gives an excellent estimate of when the bundle 
failure will take place. In our earlier work [l9| we found 
that the fiber breaking rate has a minimum at half way 
to complete collapse. However, for practical purposes en- 
ergy emission burst is a better entity to measure than the 
fiber breaking rate. 
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APPENDIX 

In the main text we have considered energy emissions 
from a slightly overloaded bundle, for which the break- 
down process proceeds slowly. Here we give a brief de- 
scription of where the energy emission minimum is lo- 
cated for larger loads. As a specific example we consider 
the model with a uniform threshold distribution. 

For a given stress a the relative number of unbroken 
fibers at the emission minimum is given by Eq.(18), 

(l-n t )nf =a i l 2 . (37) 

Since the left-hand side of (|3"Tf does not exceed 4/27, the 
maximum stress when an emission minimum is present, 
equals 

a m = 2 4/3 /9 = 0.279982. (38) 

For a given stress, a c < a < a m , Eq. (|3"T)) determines n*, 
and the corresponding value of is A/ can be calculated 
using 

nt = | + Veta4B(l.-2W*/)]> (39) 

with B = tan _1 (l/2Ve)- We have combined Eqs.(21) 
and (23) in the main text. Fig. 3 shows the result. The 
position of the emission minimum, relative to tt, de- 
creases slowly with increasing stress In this derivation 
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Figure 3: Position of the emission minimum in terms ofts/tf 
vs. excess stress for the uniform fiber strength distribution, 
computed using the continuous time approximation. The 
dashed line marks the maximum stress with an emission min- 
imum. 



the integer t is considered as a continuous variable. This 
treatment is satisfactory for slightly overloaded bundles, 
but less so when the ts is smaller. To find the smallest 
stress that produces a minimum, we consider the first 
few energy emissions. Eq. (8) for the uniform threshold 
distribution, 



E t = \Na 3 (V 3 - n^) , (40) 



and the iteration (6), no — l,n\ — 1 — a, n 2 = \_° ■, etc., 
gives 



E = a 3 (41) 




where the common factor N/6 is omitted. For all subcrit- 
ical stresses E\ > Eq, but a local minimum may occur 
already at t = 2. It is easy to show that 

E x = E 2 for a = a x = 0.278764 (45) 
E 2 = E 3 for a = a 2 = 0.273054. (46) 

For a > ax the energy emission bursts increase monoton- 
ically in size. For a 2 < a < o\ we have E 2 < E\ and 
also E 2 < E3. Hence there is a local emission minimum 
at tE = 2 in this range. 

Thus a local emission minimum is no longer present 
when the stress value exceeds <j\ . For the limiting value 
a = or, the bundle collapse occurs at tf = 8, signified by 
the first non-positive value of n t . This gives tg/tf = 0.25 
at a = 

This limiting value is a little lower than the value er m , 
the end point of the graph in Fig. 3, obtained by the 
continuous procedure above. However, at the maximum 
stress (j 1 that produces a local emission minimum, the 
value of ts/tj equals 0.230, not far from the exact value 
0.25. We conclude that the graph in Fig. 3 is very precise 
for small supercritical stresses and fairly accurate up to 
the stress a\. 



